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SOME NEW ASPECTS ON THE SUPERROTATION OF THE THERMOSPHERE 


P. W. Blum 
I. Harris 

ABSTRACT 

The motion of the thermosphere with a rotational velocity be- 
tween 10 to 20% in excess of the earth's rotational velocity has been 
deduced by King-Hele and his co-workers from the change of the in- 
clination of satellite orbits. To date no completely satisfactory ex- 
planation of the observations has been presented. In this paper we 
shall show that in the thermosphere there exists a small diurnal 
mean driving force in the eastward direction. This force has not 
previously been considered in analyses of superrotation. This 
paper presents a critical review of the observations and a theo- 
retical analysis that takes account of both equinox and solstice con- 
ditions. This work shows that the discrepancy can be resolved 
between observations and theoretical explanation in the lower 
height region where the great majority of observations were made. 
It is proposed that additional observational data are needed in the 
isothermal region for a more complete analysis. 
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SOME NEW ASPECTS ON THE .SUPERROTATION 
OF THE THERMOSPHERE 

INTRODUCTION AND REVIEW OF OBSERVATIONAL EVIDENCE 

The rotational velocity of the thermosphere as deduced from the changes of 
the inclination of the orbits of satellites during their life-time has been presented 
by King-Hele (1971) in his Figure 5. Based on 29 satellites that he has analysed, 
King-Hele's paper suggests that the ratio of the rotational velocity of the thermo- 
sphere to the rotational velocity of the earth (the superrotation ratio) increases 
linearly from 1.03 at 150 km to about 1.40 at 370 km and then decreases linearly 
from this value to a value of 0.68 at 500 km. We shall analyse the physical 
significance of King-Hele's deductions by dividing the 29 satellites into three 
groups according to height. 

The division into height groups is not arbitrary but is based on the following 
theoretical considerations: below 270 km no or very few data are available on 
the diurnal variation of the thermosphere; and, therefore, an exact theoretical 
analysis is nearly impossible, although we may estimate theoretically the range 
of the superrotation ratio (A). Between 270 and 370 km we have a farily good 
knowledge of the driving and drag forces that act in the thermosphere, and, 
therefore, we may solve the equations of motions with fairly good reliability. 
Above 370 km we would expect from theoretical considerations that no further 
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changes of velocities with height would occur due to effects of viscous drag and 
the generally assumed boundary condition that the vertical gradient of the 
velocity should go to zero. 

Table 1 shows the following for three height regions: 1) the maximum and 
minimum values of the superrotation ratio; 2) values obtained from King-Hele's 
linear fit; 3) the average values of observed A, and 4) the chi-square value for 
the test of significance of the fitted values and the average values. In the height 
region between 140 km and 270 km there are 18 satellites; thus, for this region 
we may attach greatest significance to the results of the analysis. The other 
two height regions, the second between 270 and 370 km and the third above 370 
km, contain only four and three satellites respectively; therefore, obviously any 
conclusions will have less significance. We suggest that King-Hele's interpre- 
tation of his analysis is not the only possible one. The following interpretation 
is also equally valid: 

1. There is reliable evidence that in the height region between 140 km and 
270 km the thermosphere rotates faster than the earth. The value of 
the supperotation ratio is probably between 1. 1 and 1. 3. The mean 
value deduced from observations is 1. 17, but when the observations 
are weighted according to their observational errors, the mean is re- 
duced to 1. 14 — resulting mainly from the small error of Cosmos 
1969-94b, which yielded a superrotation ratio of only 1. 06. If it is 
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Table 1 


Height 

Region 

km 

No. of 

Observed King-Hele 

Average 

Chi-square Test 

Satellites 

Max 

Min 

Fit 

King-Hele 

Average 

150-270 

18 

1.4 

1.0 

1. 05-1.23 

1. 15 

1. 76 

1.84 

270-370 

4 

1.6 

1.1 

1. 25-1.40 

1.26 , 

2. 03 

2.85 

370-500 

3 

1.35 

0.71 

1.4 -0..69 

1.04 

— 

— 


assumed that the superrotation rate is not height dependent but has a 
constant value in this height region, a chi-square significance test 
shows that such an assumption does not result in any significant change 
in the probability of goodness of fit as compared with King-Hele's fitted 
linear height dependence. 

2. In the height region between 270 and 370 km King-Hele has analysed 4 
satellites. He deduces for this height region an increase of the super- 
rotation ratio from about 1.2 at 270 km to 1.4 at 370 km. Again we 
must conclude that in this region the superrotation ratio A is larger 
than unity, probably between 1. 1 and 1. 4. Due to the very limited 
number of satellites analysed, a deduced height dependence of the 
superrotation ratio A becomes questionable from a statistical point of 
view, although the observations do not exclude a height dependence. 

3. In the third height range above 370 km King-Hele deduces decreasing 
superrotation ratios. At 500 km he suggests a ratio of 0. 68. Only 
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3 satellites were analysed in this height region. The arguments for an 
assumed height dependence — considering the statistics — are even 
weaker than in the second height region. 

It is very difficult to reconcile a height dependence of the superrotation ratio 
above 370 km with theoretical considerations. On the other hand, we shall show 
that there latitudinal dependence of the superrotation ratio is possible by theo- 
retical considerations with decreasing superrotation ratios for satellites having 
large inclinations. In this context it should be observed that 1967-42A (Ariel 3) 
resulted in a superrotation ratio of 0. 7 at 500 km. This satellite has an inclin- 
ation of about 80° . For these reasons we suggest that the data available so far 
do not justify the deduction that the superrotation ratio decreases so markedly 
above 370 km as suggested by King-Hele. 

Figure 1 shows the observations, King-Hele' s fit and our suggested mean 
values for the three height regions. 

THE DIURNAL AVERAGE AZIMUTHAL DRIVING FORCE 

The equations of horizontal motion in the thermosphere show that a diurnal 
mean motion of the thermosphere can result from four effects: 1) A diurnal 
mean-zonal driving force; 2) The variation of the drag forces, especially ion 
drag with local time out of phase with the zonal velocities; 3) Ion motion driven 
by meridional on radial electric fields; 4) Mean meridional motions that inter- 
act through the Coriolis force on the zonal motions. 
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Rishbeth has recently (1972) reviewed the various possibilities of exjriaining 
theoretically the observed superrotation and concludes that more data are re- 
quired for a well-founded theoretical explanation. 

In all previous treatments it has been assumed that the mean diurnal zonal 

i * . . ' * • * .r 

driving force vanishes, as this is a consequence of the Jacchia model (1971). It 
can be shown that based on the Jacchia density distribution and the observation 
of thermospheric temperatures by incoherent radar back scatter measurements 
a small eastward diurnal average driving force results. 


The zonal driving force is given by (Blum and Harris, 1973) 

e - R / x 9, 9T T' 3M \ 
d sin0 Mtoj ( T 9r " P 3r . M* d r ) . 

Where p is the pressure, T the temperature, P the density, M the mean molec- 
ular weight, t the local time, 0 the colatitude and R the universal gas constant. 
We shall neglect changes of the mean molecular weight with local time. Then 
the expression (1) reduces to 


R /3 dT \ 

f = — T — In p + — ) . • (2) 

d sin0 Mrcu \ dr dr ) 

We are interested in the diurnal mean (denoted by < > ) of the driving force. To 
this mean the term 3 T/3r does not contribute, therefore 


<f d > = - 


R ^ 9 

< T — In p > 

sina Mur. 3r 


(3) 
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Writing T and log p in Fourier components 

m 

T = ^ a k cos a) k (r - 
k=0 

In p = 2 b k coscok(T-p^b (4) 

where p (J) and p(£) are the phases of the Fourier components of the temperature 
and log density. We obtain for the diurnal mean, 

n rn i 

<f d> = SrfSTS g - V k “ sin (p<P - P<{>) (5) 

Generally only the term with k = 0 has been considered and no mean diurnal zonal 
driving forces results. We shall estimate the contribution of higher Fourier 
terms by using the observationally verified phase difference of about 1. 5 hours 
between density and temperature in the thermosphere. For this simple estimate 
we shall assume that the phase difference of 1. 5 hours between temperature and 
density arises only from the first, or diurnal, Fourier components. 

C f d > = - * aj bj sin (22.5°) (6) 

u sm0 Mr 1 1 

Substituting the numerical values for the diurnal amplitudes of temperature 
and lnp found in the thermosphere, we obtain for the mean driving force at 
300 km 

f d/o> = 14 - 5 m / sec 

We have calculated the value of the mean zonal driving force at a height of 300 
km accurately from the Jacchia model that was modified by assuming that the 
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temperature peaks about 1. 5 hours later than the densities. This result is shown 
in Figure 2. It is seen that for this height the mean diurnal zonal force divided 
by co, the earth's angular velocity, is at the equator about 7 m/sec. At a height 
of 500 km this force has increased by a factor of approximately 3. While this 
force induces some mean eastward velocity, it alone is insufficient to explain 
King-Hele's observations. 

In the lower thermosphere we have no data on the phase difference between 
density and temperature although a much larger phase difference than 1. 5 hours 
is probable. The maximum eastward driving force would result if the temper- 
ature peaked 6 hours after the density. This would increase the mean diurnal 
zonal force by a factor of 3 as compared with the forces used in our computation. 
As the drag forces in the lower thermosphere are small, the increase of the 
superrotation factor due to the possible increases of the driving forces below 
250 km may be considerable. As no data on the phase difference below 250 km 
are available, we have also used for the lower height region a 1. 5 hour phase 
difference. 

THE EFFECT OF DAY-TO-NIGHT VARIATION OF ION DENSITIES ON THE 
MEAN ZONAL VELOCITY 

The particularly simple case of equinox conditions at the equator lends itself 
to simple estimates of the possible influence of the ion density variation on the 
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mean zonal velocities. The equation of horizontal zonal motion, V , is 

dV<^ 3 2 v^ * 

— +2cjcos 0V (6) - V dz 2 + d ion V ( ^> - f d 


(7) 


where 0 is the colatitude, r\ the kinematic viscosity, w the velocity of the earth, 
dj on the ion drag coefficient, and f d the zonal driving force (Blum and Harris, 
1973). 


In the height range where ion drag dominates over viscous drag, the equation 
for the diurnal mean motion at the equator becomes 

<d i0 n v( ^> - <f d > < 8 > 

or explicitly in terms of Fourier coefficients up to the diurnal terms 

v(*> - 4-£iV rt cos(p (d)_ pl V) ) (9) 

d'*» 2 d Q V l H l 

where d 0 is the mean diurnal ion drag, dj, V, the diurnal amplitudes of ion drag 
and zonal velocity respectively, p^ and p'^ the phases of the diurnal components of 
ion drag and zonal velocities respectively. 


Based on the results of Blum and Harris (1973) for the numerical values of 
the ion drag coefficient and the zonal diurnal motion, we may estimate the re- 
sulting mean diurnal zonal motion. We shall use the following numerical values 
in Table 2 as representative. We obtain the tabulated results for the mean zonal 
motion. Obviously the values of A are too low for an explanation of the 
observations. 
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Table 2 


d o 

d l 

cos (pj^ - p^ v) ) 
fo 

y(V’) 

yW 

o 

A 


220 km 


1.4 


1.4 
~- 0. 5 

2. 2 m/sec 
115 m/sec 
30. 3 m/sec 
1.07 


340 km 
9.6 
2.27 
0.5 

10.1 m/sec 
56 m/sec 
4.4 m/sec 
1.01 


Rishbeth (1971) has suggested that due to polarization effects in the F region 
the ion drag is reduced during the day by a factor R day of about 0.8 and by night 
by a factor R ni of 0. 2. While Rishbeth's suggestion is not completely verified 
or theoretically founded, (Volland, 1971), we may easily estimate its effect on 
the mean diurnal motion for the above simple case. We shall have to assume that 
even with Rishbeth's modification the dominance of ion drag over viscous drag 
is maintained. This assumption may yield an overestimate of the effect. 

Denoting the ion drag coefficient change by R day for the daytime values and 
R ni for the night values, it is seen that Rishbeth’s modification amounts to re- 
vised ion drag coefficient (d 0 * d ; *) as follows: 

d o* *2^° ^day + + ~ d l ^day ” 

d l * - “ d o ^day ” ^-ni^ + *2 d l ^day + ^ni) (10) 

with the phases of d and d l equal. 
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We shall assume that the modified diurnal amplitudes of the zonal velocity 
vary inversely with the mean ion drag. Thus, we obtain for the ratio of the 
modified mean zonal velocity V q ( r ) to the unmodified mean zonal velocity V Q the 
expression 


/-“S 

oc 

'w' 

o 

> 

2 

1 

2 

j + 4 d o R day R ni 

V o 

R day + Ki 

1 + d l . 2 R day “ R ni 
d o n R day + R ni 


* d l R day + R ni 


It is seen that the ratio of modified to unmodified mean zonal velocities at 
the equator is not very sensitive to the particular values of H day and R ni but only 
to their difference. This difference has, according to Rishbeth’s suggestion, a 
value of 0. 6. This yields the estimate for the ratio of the revised mean zonal 
velocity to the unmodified mean zonal velocity (according to equation (11)) be- 
tween 4 and 5. Thus, from this simple estimate, we obtain mean zonal velocities 
that are in the range of 20 to 30 m/sec in the isothermal' region. Rishbeth's 
modification is not directly applicable to the lower thermosphere. 

Rishbeth, who has also calculated the effects of his modification on the mean 
zonal velocities, obtained a factor of 25 instead of our value of 4 to 5. Thus, we 
cannot confirm his estimates. We shall present later the results of an accurate 
calculation of the mean zonal velocities with Rishbeth’s modification. This exact 
calculation yields a ratio of the modified to the unmodified cases of about 6. As 
our unmodified velocities are higher than the previous estimated values of 
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4.4 m/sec (they are 10 m/sec), we obtain mean zonal velocities at the equator 
of about 55 m/sec as shown in our Figure 3. 

DEFINITION OF THE SUPERROTATION RATIO 

From a given distribution of the global zonal velocities we may determine 
the superrotation ratio A that is observed from the change of the inclination, i, 
of the satellite as analysed by King-Hele. 

As the zonal force acting on the satellite is proportional to p V^\, we have 
to calculate as a first step the weighted mean zonal motion V 1 ^ defined by 

V(*) (0,Z) = / V^)(t,0,Z)p(t,0,Z)dt 
where the integration is over one day. 

The weighted mean motion is generally somewhat less than the mean diurnal 
velocity V 0 (ip) . In the following we shall use V^) instead of vfa) as the difference 
is not very considerable. 

The observed change of the satellites inclination is accumulated at all lati- 
tudes during the satellites motion. In line with King-Hele’ s definition of the 
superrotation ratio A we calculate 
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where v ea is the earth's velocity at the equator, 6 the colatitude, i the inclination 
of the satellite's orbit (identical with the maximum latitude reached by the 
satellite). 

The integration is extended over all latitudes from -i to +i. The integral in 
the demoninator may be evaluated in a closed analytic form. 

DISCUSSION OF RESULTS 

The method of integration of the horizontal equations of motion of the thermo- 
spheric wind field has been described by Blum and Harris (1973). This method 
includes all the non-linear terms of the equations. The equations were integrated 
with Rishbeth's suggested modification of the ion drag coefficient due to F region 
polarization fields. The modification that was originally only made for the 
equator was generalized to all latitudes. The generalization assures continuity 
everywhere and the ion drag at the poles remaining unmodified. 

Figures 3 and 4 show the results of the computations of the mean zonal 
velocities for equinox and solstice conditions respectively for a height of 300 km 
for Rishbeth's suggested values of the reduction of the ion drag by a factor of 
0. 8 (Rday) by day and a factor of 0. 2 (R ni ) by night. Shown also are the non- 
linear solutions for the unmodified ion drag. No convergent non-linear solution 
was obtained for solstice conditions with Rishbeth's suggested modification; so 
only the linear solution is given. It is seen that the mean zonal velocities de- 
crease faster than the earth's rotational velocity with increasing latitudes. An 
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integration over latitudes, as described above, has been performed in order to 
obtain the superrotation ratio as a function of the inclination of satellite orbits. 
Thus, although relatively large mean zonal velocities are obtained at the equator 
(about 112 m/sec at solstice and 57 m/sec at equinox) with Rishbeth's modifi- 
cation, the effective values of the superrotation ratio so deduced are much 
smaller than the mean equatorial velocities except for very low inclination 
satellites. 

The superrotation ratios are illustrated in Figures 5 through 8 as a function 
of the maximum latitude that is reached by the satellite, i. e. , its inclination. 

It is seen that even Rishbeth's modification results in a decrease of the superrota 
tion ratio with increasing latitudes. Around a narrow band of latitudes near the 
equator the superrotation ratio at 340 km is 1.13 for both solstice and equinox 
conditions. For higher inclination satellites it decreases to 1. 09. For solstice 
conditions at 240 km a value of A = 1. 07 is obtained for the equatorial zones; 
it decreases for higher latitudes to 1. 05. At equinox at a height of 240 km 
Rishbeth's modification yields A = 1. 13 for the equator, and it decreases to 1. 05 
as the latitude increases. The rapid decrease with latitude for the non-linear 
solution is partly due to the increase of the mean equatorward meridional wind 
(Blum and Harris, 1973) that causes through the Coriolis term a westward 
driving force. 
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Radial and meridional electric fields could cause an ion motion in the zonal 
direction, thus effectively changing the ion drag. They may even cause the 
neutrals to be dragged by the ions; i. e. , the ions could have a larger velocity 
than the neutrals. If such fields are included in the equations of motion, they 
would be independent of the neutral velocity and could be included in the driving 
force on the right hand side of the equations of motion. Arbitrary assumptions 
of the direction, the magnitude and the time dependence of these fields could ex- 
plain any mean zonal wind field. It seems that no evidence exists regarding the 
existance of electric fields that could explain the high values of the superrotation 
ratio found deduced by King-Hele for the heights near 370 km. In fact, Harper 
(1971) has deduced meridional electric fields at Aricebo that would cause a west- 
ward force acting on the thermosphere. 

CONCLUSIONS 

1. The superrotational velocity of the thermosphere below 270 km that has 
been determined by King-Hele from satellite drag analysis yields a 
ratio of about 1. 15 to the earth's rotational velocity. The height de- 
pendence below 270 km suggested by King-Hele can neither be confirmed 
nor negated by statistical tests of significance. 

2. Between 270 and 370 km there also exists an eastward rotation of the 
thermosphere. Satellite drag analysis seems to indicate a superrotation 
ratio of 1. 25 with a maximum of 1. 4 at 370 km. Statistically no definite 
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determination of the superrotation ratio— and even more so its height 
dependence— can be made as not enough data are available. 

3. Above 370 km only three satellites have been analysed. The results 
were interpreted by King-Hele as showing a decrease of the superro- 
tation ratio with height to a ratio less than unity at 500 km. It is sug- 
gested that the deduced decrease may be a latitudinal effect rather than 
a height effect. 

4. A theoretical treatment involving the integration of the equations of 
motion for the thermosphere yields the following results: 

a. Below 270 km superrotation ratios that are compatible with the 
observations may be deduced. 

b. Without extremely hypothetical assumptions like strong radial 
electric fields no superrotation ratios in excess of 1. 17 can be 
calculated. This makes it impossible to explain theoretically 
King-Hele's interpretation of the observations between 270 and 370 
km. 

c. A theoretical treatment will not result in a decrease with height of 
the superrotation ratio above 370 km unless the theoretical assump- 
tions regarding the transition region at the exobase are abandoned. 
This absence of a height dependence of the superrotation ratio above 
370 km is not in conflict with the calculated mean driving forces. 
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These are due to the phase difference between density and temper- 
ature in the thermosphere and contribute to the observed 
superrotation. 

d. The polarization effects in the F region suggested by Rishbeth in- 
crease the superrotation ratio by a factor 3 to 6, depending upon 
latitude and altitude. This is significantly less than Rishbeth’ s 
determination of the factor of 25 at the equator. This factor is 
about 6 at the equator but only about 3 for satellites having high 
inclinations. 

e. In the height region between 340 and 500 km there results an in- 
significant height dependence of the superrotation ratio. 

f. Generally the superrotation ratio is larger at solstice than at 
equinox. Using Rishbeth's suggested values of R day = 0. 8 and 
R ni = 0. 2 there results at solstice a superrotation ratio of 1. 12 
at 300 km. Above 400 km this ratio is increased to 1. 15 while in 
the lower height region at 240 km it is about 1. 09. 

g. The effects of the non-linear terms of the equations of motion re- 
duce the superrotation ratio, thus demonstrating that it is not a 
non-linear effect. 

The extension of the modification of the ion drag coefficient as suggested by 
Rishbeth to mid and high latitudes does not yield a superrotation ratio in excess 
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of 1. 15. We can theoretically obtain superrotation ratios between 1. 1 and 1. 15 
by taking into account the mean eastward driving forces, the modification of the 
ion drag due to polarization fields in the F region and the motions at solstice 
conditions (which generally have higher mean eastward velocities than at equinox) . 
Higher superrotation ratios cannot be explained without additional assumptions 
like radial or meridional electric fields having a particular time and space dis- 
tribution. Such fields are hypothetical. 
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FIGURE CAPTIONS 


Figure 1. King-Hele's results of the height dependence of the superrotation 
ratio based on 25 satellites. The fully drawn lines show the re- 
sults of an assumed constant superrotation ratio separately 

for each of the three height regions defined in the test. The sta- 
tistical reliability of such distribution of the superrotation ratios is 
nearly equal to that of King-Hele's suggested height dependence. The 
analysis gives to each observation a weight according to the ob- 
servational error given by King-Hele. 

Figure 2. The mean azimuthal force for summer solstice conditions as a 

function of latitude, for the Jacchia model with a phase difference 
of 1. 5 hours between temperature and density, at the altitude of 
300 km. 

Figure 3. The mean zonal velocity as a function of latitude at the height of 

300 km for equinox conditions. ( ) is the linear solution with 

Rishbeth's modification reducing the ion drag by day by the factor 

R day - 0. 8, and by night with a factor R nj = 0. 2. ( ) is the 

non-linear solution with Rishbeth's modification, and ( ) is 

the non-linear solution without Rishbeth's modification. 
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Figure 4. 
Figure 5. 


Figure 6. 
Figure 7. 
Figure 8, 


Same as Figure 3 but for solstice conditions and no non-linear 
solution for Rishbeth' s modification is given. 

Superrotation ratio as a function of inclination of a satellite’s orbit 
for equinox conditions at the height of 240 km. ( ) is the re- 

sult for linear solution with Rishbeth' s modification, 

the non-linear solution with Rishbeth’ s modification, ( ) is the 

linear solution without Rishbeth' s modification, (....) is the non- 
linear solution without Rishbeth' s modification. 

Same as Figure 5 but for a height of 340 km. 

Same as Figure 5 but for solstice conditions. 

Same as Figure 6 but for solstice conditions. 
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